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1. Introduction 

In 1978, Lozi introduced in a short note pL] a two-dimensional map the equations and 
attractors of which resemble those of the celebrated Henon map [2] . Simply, a quadratic 
term in the latter is replaced with a piecewise linear contribution in the former. This 
allows one to rigorously prove the chaotic character of some attractors [3J and a detailed 
analysis of their basins of attraction [4J. The equations for the iterated Lozi map L{x, y) 
read 



where a, b are real non-vanishing parameters. Inside the region where the orbits remain 
bounded, the Lozi map may present both regular and chaotic behavior. 

The purpose of this article is to show that the Lozi map also exhibits a special 
type of bifurcation recently found in a class of one dimensional discontinuous piecewise 
linear maps [5]. A similar phenomenon in continuous maps was also described in [6]. 
The main feature of this bifurcation is the presence of a continuum of neutrally stable 
cycles. Thus, at the very bifurcation point, the attractor in phase space contains an 
infinite set of regular orbits. As a function of the parameters the transit may take place 
from regular to regular as well as to chaotic regimes. 

In maps in one dimension this kind of transitions happens whenever a piece of an 
iterate of the map becomes co-linear with the identity function, the bisectrix, in the 
cobweb for a particular value of the parameter. We will extend this explanation to 
maps in two dimensions and discuss its connection with the theory of border collisions, 
which is commonly used to analyse these kinds of maps [3, El E]- Very recently, this 
question has been considered from a more abstract point of view [10] and the name 
degenerate bifurcation has been used. As this term may sound polysemous we refer to 
this phenomenon as bisecting bifurcation in the rest of the paper, which seems more in 
accordance with the geometrical interpretation advanced. 

An overview of the results on the Lozi map necessary for later reference is presented 
in section [2j Bisecting bifurcations are introduced and analysed in section [3j They are 
then given a geometrical explanation in section |4j In section [5] we further illustrate how 
the dynamics of the map changes with parameter a and establish a connection with 
the theory of border collisions. Section |6] contains some final comments on somewhat 
related topics. 

2. Dynamics of the Lozi map: isolated attractors and stability 

It is customary to start the study of a dynamical system with a catalogue of its fixed 
points, periodic attractors and so on, paying special attention to their stability and 
evolution with the system parameters. We start by collecting some of these features for 
the Lozi map which will be referred to in the following sections. 



Xn+l = l+Vn- a,\x. 
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For stability considerations we need the Jacobian matrix of L{x, y) which reads 

/(^,y)=(-*' j). (2) 

Notice that J{x, y) depends on the point of the orbit only through sign(a;). Accordingly, 
we denote its values as J+ = J{x > 0,y) and J_ = J{x < 0,y). Furthermore, since 
det J± — —b, we will only consider maps with \b\ < 1, in order to have non-expanding 
systems. In particular, the maps with b — ±1 are area-preserving. 

It is a simple exercise to see that for |6| < 1 the Lozi map has no fixed points if 
a < b — 1. It has Pi as the unique fixed point if 6 — 1 < a < 1 — b, and gains an additional 
fixed point, P2, ii a > 1 — b, with 

Vl + a — 1 + a — b/ VI — a — o 1 — a — b/ 

Furthermore J{Pi) = J+, J{P2) = J-- Only Pi may be stable and this 
happens for systems with parameter values (6, a) inside the triangle with vertices 
(1, 0), (—1, 2), (—1, —2) in parameter space. 

As for isolated cycles of period T = 2, it is easy to see that they exist, if we keep 
|6| < 1, for a > 1 — 6. Their elements are {zi,bz2) and {z2,bzi), with 
1 — b + a 1 — b — a 

" (1 - 6)2 + a2 " (l-6)2 + a2- 

They are stable for parameter values (6, a) inside the triangle with vertices 

(0, 1), (1, 0), (1, 2) in parameter space. 

The analysis so far is standard and may be pursued for more complicated attractors. 
Wc have recalled the previous facts just to better understand the bifurcations which will 
be studied in the next section. 



3. Bisecting bifurcations and attracting sets 

The fixed points and cycles mentioned above were isolated points in the phase space. 
However, for particular values of the parameters a and b we can also find continuum 
sets of periodic points which act as attractors for certain regions of phase space. We 
commence by explaining the approach we follow to locate algebraically which values of 
a and b give rise to such infinite sets of simultaneous cycles. We solve first the cases of 
periods two and four which we deem illustrate sufficiently the procedure. After that, 
we completely determine the cycles' elements. Period three and five are then analyzed 
in the same way and shown to produce a different pattern. The stability of the cycles 
is also resolved. Higher order periodic attractors may be studied along the same line, 
of course at the price of more involved algebra. 
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0, rank I « ^ ' 7 , i 1 = 1' (») 
1 — a sign(x2) 



3.1. Continuum of period— 2 attractors 

Let {(xi, (x2, 1/2)} be a period— 2 cycle of the Lozi map. Next, instead of ([T]) we use 
the equivalent second order difference equation 

Xn+l = 1 ~ O-IXnl + bXn-l, (5) 

with the obvious corresponding change in the initialization of the iteration. The 
periodicity condition for this system conveys 

X2 = 1 — + bx2, (6) 
xi = 1 — a|x2| + bxi. (7) 

The corresponding 2— cycles of the Lozi map will be {(xi, 6x2), (x2, bxi)}. 

When the system of equations above is compatible but indeterminate then there 
exists an infinity of solutions and hence of period— 2 cycles. From an algebraic point of 
view this means that the coefficient and the augmented matrices must both have rank 
one: 

a sign(xi) 1 — 6 
1 — 6 a sign(x2) 

which yield the constraints: a > and either 6 + a = lor6 — a = l. Correspondingly 
there are two one-parameter families of 2— cycles in M^, {^i, S2} and {Ui, U2} with 

= (x, 6(1 — ax)/a), ^2 = ((1 — ax)/a, 6x), (0 < x < 1/a, 6 = 1 — a) (9) 
f/i = (x, -6(1 + ax)/a), U2 = {-{I + ax)/a,bx), (-1/a < x < 0, 6 = 1 + a). (10) 

For the S"— family we have J{Si)J{S2) = whose eigenvalues are (a^ + 26 ± 
|a| + 46)/2. In the critical case, a + 6 = 1, they reduce to the values 1 and (1 — a)^, 
and hence these orbits are neutrally stable whenever < a < 2. Then the continuum 
they form is an attractor in phase space. 

For the [/—family we have J{Ui)J{U2) = Jl, whose eigenvalues in the critical case, 
6 — a = 1, are 1 and (1 + a)^ > 1, with a > 0. Therefore these orbits are unstable. 

3.2. Continuum of period— A attractors 

The search of sets of period— 4 orbits yields the following results. First, the constraints 
on the parameters are a = 1 + 6, and 0<6<1, l<a<2, for the rank of the 
corresponding matrices to be equal to three, allowing in this way the existence of a 
one-parameter family of solutions. Second, if (xj,?/j) denotes the cycle elements, with 
i = 1, . . . , 4, then 

xi = /i + 26/3, X2 = /i — (6 — l)/3, X3 = — /i, X4 = — /i + (6 + (11) 
yi = bx4, y2 = 6x1, ys = 6x2, 1/4 = bx^, (12) 
where we have used the real parameter < < 1 and the definition /3 = 1/(1 + 6^). 



One can readily deduce from the expression of the cycle elements (11) that the 



X— components are alternatively positive and negative. Therefore the corresponding 
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Jacobian matrix reads simply [J^J^Y, whose eigenvalues are 1 and (1 — a)^ in the 
critical case. Hence, these period— 4 orbits are neutrally stable and their union acts also 
as an attractor. 

3.3. Period— 3 attractors 

Following the same algebraic procedure as in the preceding subsections, we have 
determined that an infinite set of cycles of three elements exists only if a = 1 and b = —1. 
However, in this case the rank of the corresponding coefficient and augmented matrices 
is equal to one and, as a consequence, we get a continuum of solutions depending on 
two parameters. 

A phase portrait which includes the period— 3 orbits is shown in Figure [3j Orbits 
of period— 3 fill completely the white triangle at the center of Figure |3j with right angle 
vertex at (0,0). The acute angle vertices are located at (1,0) and (0, 1). The elements 
of the cycle are given by (^, —rj), (1 — ^ — 77, — ^) and (77, — 1 + ^ + rj), in terms of the 
parameters < ^,,1] < 1. The remaining white areas correspond to regular orbits of 
higher periods. Finally, the black and grey areas stand for chaotic trajectories. 

The two eigenvalues of the Jacobian matrix of L^^^{x,y) evaluated on the orbit 
elements are equal to one, which stands for double neutral stability. Hence, the set of 
orbits of period— 3 is two-dimensional in phase space, at variance with period— 2 and 
—4, and is embedded into the sea of chaotic trajectories. 

3.4. Period— 5 attractors 

The situation with period— 5 attractors is similar to the one with period— 3. We have 
found two two-parameter sets of orbits with neutral stability. They emerge when b = —1 
and either a = {1 + V^)/2 = (the golden ratio) or a = (1 — -\/5)/2. Table [I] 
collects the elements {xj,?/,}, with i = 1,...,5, of the orbit. We use the notation 
(f = (y/b — l)/2 = — 1 = 1/0. The two parameters u and v which define the family of 
period— 5 orbits vary in the pentagon of vertices {(v^^, 0), (y?, 0), (y?^, ip^), (0, ip), (0, ip"^)} 
in parameter space, provided a = (p. Else, if a = —ip, the two parameters vary in the 
pentagon of vertices {(0, 0), (1, 0), (0, 1), (1, 0), (0, 1)}. As far as stability is concerned, a 
direct calculation shows that the Jacobian matrix of L^^\x,y) for the parameter values 
just quoted is indeed the unit matrix which means double neutral stability. 

3.5. Bifurcation diagrams 

It is customary to summarize the behavior of a map with varying parameters by means 
of bifurcation diagrams. In principle there are two somewhat extreme ways to construct 
these diagrams. The ideal mode to proceed is clear: one determines by a theoretical 
analysis the periodic points and their stability and plot them as functions of a parameter. 
This is often impossible to carry out completely and then one has to resort to a more 
experimental way: iterate the map for a number of initial conditions, discard a transient 
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in the trajectories and plot then a component as a function of the parameter. The 
unavoidable limitations of the last approach are clear and one has to take special care 
to identify the results obtained. In particular, a densely filled bar in a diagram is usually 
taken cLS cL SI mature of chaos. 

In the Lozi map, when enough initial conditions are considered for the numerical 
simulation, the continua of 2— and 4— cycles discussed above may act as an attractor 
and show up also as a vertical line in the bifurcation diagram. This is clearly seen in 
Figure [T]for fixed 6 = 0.1. To construct this diagram especial care has been payed to 
sweep a large number of initial conditions. For every trajectory a large enough transient 
is discarded. Also, if a component increases beyond a given threshold the trajectory is 
considered as unbounded. Figure [T] exhibits a bifurcation from period— 1 to period— 2 at 
a = 0.9 and another one from period— 2 to chaos at a = 1.1. This numerically obtained 
diagram illustrates clearly that these two bifurcations are mediated by the phenomena 
described in sections |3 . 1 1 and |3l2] respectively. A further feature of Figure [T] is that in the 
approximate intervals 1.1 < a < 1.237 and 1.237 < a < 1.406, the trajectories jumps 
between the chaotic bands in a cyclic manner, with periods 4 and 2 respectively. 

Similarly, in Figure |2| here for fixed a = 1.5, two bifurcations occur at 6 = ±0.5, 
again mediated by the same mechanism. The system is chaotic for b = 0.1 and a > 1.1 
in Figure [T| and a = 1.5, —0.5 < b < 0.5 in Figure |2} The presence of these vertical 
segments corresponding to the described attractors composed of infinite sets of periodic 
orbits is the most salient feature of the phenomenon we deal with here. We term it 
bisecting bifurcations for reasons which will appear clear in the next section. Eventually, 
cyclic jumping between the chaotic bands of the diagram takes place with period 2 and 
4 to the left and right of the figure, respectively. 

We gather the information discussed so far in Figure |4] where the behavior of the 
Lozi map in different regions of the parameter space is indicated. Trajectories for |6| > 1 
are unbounded. Numerical simulations indicate this is also the case for the northern and 
southern parts of the diagram. The identification of the chaotic area has been achieved 
by computing numerically the Lyapunov exponents [H], thus the upper border of chaos 
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Table 1. Elements of the period— 5 orbits {xi^yi}, {i = 1,...,5). The upper block 
entry is for a — (j). The real parameters m, v take values inside and on the pentagon 
of vertices {((ys^, 0), {ip, 0), (y;^, (/?^), (0, (p), (0, ip^)} of parameter space. The lower block 
entry is for a ~ —Lp and the parameters are inside and on the pentagon of vertices 
{(0, 0), (1, 0), (0, 1), (1, 0), (0, 1)} in parameter space. In both cases it is 6 = —1. All 
the orbits elements stay in the first quadrant of phase plane. 
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is approximate. The stable fixed points of L and L^^l, are included in the graphics. The 
borders of their stability region will be analysed in detail later on. The diagram by no 
means exhausts the variety of behaviors the trajectories of the Lozi map may develop 

m. 



4. A geometrical explanation for the bisecting bifurcations in the Lozi map 

In [5j, a geometrical explanation for the bisecting bifurcations in piecewise one- 
dimensional maps is provided. Here we develop a generalization for the two-dimensional 
Lozi map, although the graphical illustration is more involved. Figure |5] represents the 
surface f^'^^{x,y) and allows us to appreciate how the simple shape of f{x,y) in ([l]) 
becomes convoluted after iteration. 

Notice that for the Lozi map, the existence of a bisecting bifurcation in the variable 
X implies its presence in the variable y, and viceversa. This is due to the fact that y 
is just a one iteration delayed and re-scaled version of x. This observation will render 
easier the design and interpretation of some three dimensional diagrams. 

From an algebraic point of view, we will find an infinity of period-ra limit cycles if 

/["l(x,y)=x, g^'^\x,y)=y, (13) 

for some range of (x, y) values, as illustrated for n = 2 to 5 in section [2] Clearly, this 
is the condition of fixed point extended to a continuous set of points in phase space. 
Furthermore, this is only a necessary condition for the observation of the infinity of limit 
cycles since its stability is not assured by it. 

In one-dimensional maps it is a common practice to construct and study trajectories 
using the cobweb diagram. In a two-dimensional generalization we can visualize /t"' (x, y) 
and g^"^{x,y) as surfaces to which a point p = {xm,ym) of the phase plane is projected 
upward in order to find the new Xm+i = f^'^\xm-,ym) and ym+i = g^"'K^m,ym) on the 
vertical axis. This procedure is illustrated in Figure |6] for the particular case n = 1. 
The horizontal plane stands there for the phase plane. The vertical axis is common for 
both surfaces /^"^(x, y) and g^"'^{x, y), which have been sketched only partially. Using the 
cobweb diagram technique, the values Xm+i and ym+i are sent back to the phase plane. 
Thus, the value Xm+i is first projected toward the bisectrix hne /'""'(x, 0) = x and then 
down to the x-axis. Similarly, ym+i is first projected toward the bisectrix g^'^^{0,y) = y 
and then down to the y-axis. This gives the point q = (xm+i,l/m+i) on the phase plane 
of Figure [6j 



Let us now interpret condition (13) in geometrical terms. The presence in a 
two-dimensional map of a bifurcation mediated by an infinite set of neutrally stable 
cycles implies the existence of a segment, or set of segments, in the phase plane whose 
projections upward on the surfaces /^"''(x, y) and (7'"'(x, y) are contained in the bisecting 
planes f^'^^{x,y) = x and g^"'^{x,y) = y respectively. This explanation justifies the name 
bisecting bifurcation we have given to the phenomenon we are discussing. 



In the case of the Lozi map, / and g allow explicit solutions for condition (13) to be 



found, at least for low values of n. Next we focus on the particular case of the bifurcation 
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occurring at a = 1 — b with a > 0. In this case the hmit cycles at the bifurcation point 
are all period— 2 cycles, i.e. fixed points of L'^^ {x, y). In Figure[7]we have represented the 
surface /'^^(x, y) and the generalized cobweb procedure described above. For the sake of 
clarity, any analogous representation concerning g^'^^{x,y) has been omitted. The plot 
has been generated for b = 0.1 and a = 0.9. Any initial condition (xo,?/o) lands, after a 
certain transient, on the segment y = b{l — ax) /a with < x < 1/a (red line in Figure 
I?]). Every point of this segment is a fixed point of L^'^^{x,y). When projected to the 
surface /^^'(x, y) they graze the bisector plane /'^^(x, y) = x and then come back to their 
position in the phase plane. Since this segment is invariant under the transformation 
L^'^^(x,y), any point in the segment is an element of a period— 2 cycle of L{x,y). 

In summary, from a three-dimensional geometrical perspective, we will find an 
infinite set of limit cycles of period n from the Lozi map when the projections on the 
phase plane of the two curves produced by the intersections of i) the surface f^^^{x,y) 
and the bisector plane z = x, and ii) the surface g^'^^{x, y) and the bisector plane z = y, 
overlap in a certain range of the variables x,y. Whenever the intersection is not a linear 
segment but an area, the system has double neutrally stable orbits, which is the case of 
the period— 3 and —5 orbits we have studied. 

5. Bifurcations and phase space 

In this section we discuss the appearance of the bisecting bifurcations in connection with 
the changes in the stability of other attractors. This will further clarify the meaning 
of the vertical segments in bifurcation diagrams like the ones in Figures [T] and [2j To 
be specific we fix 6 = 0.5 and plot in Figure |8] twenty regular attractors in phase space 
corresponding to parameter values a = 0.1 n, n = 1, 2, . . . , 20. 

The diagram allows us to appreciate how the fixed point attractor evolves with 
increasing a for n = 1, . . . ,4. At a = 1 — b = 0.5, i.e. n = 5, the first bisecting 
bifurcation takes place. Observe that this is precisely where the neutrally stable fixed 
point Pi (a, 6) in (|3| changes to be unstable. We note in passing that this type of 
behavior sounds familiar from the analysis of elementary systems. For example, in the 
most conventional period-doubling cascade that would mean the birth of a period— 2 
attractor. Here something of this sort is also observed because 2— cycles exist to the 
right of a = 0.5. The difference is that the transition is mediated by a continuum of 
neutrally stable cycles of period— 2 which fill up the segment in the figure. 

When we continue to increase a the system is periodic with stable period— 2 cycles. 
These have been represented in the figure for nine values of a, namely n = 6, . . . , 14. 
Then, at a = 1 + 6 = 1.5, the second bisecting bifurcation takes place. Again these 
parameters values punctuate the transition from stable to unstable for the 2-cycles 
mentioned in section |2} The two segments of the attractor are now built up by neutrally 
stable cycles of period— 4, and they are apparent in Figure |8| Values a > 1.5 convey 
chaos (not shown) with unstable cycles in it (n = 16, . . . , 20). 

The case considered illustrates, for fixed b and varying a, how the Lozi map follows 
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a route to chaos in three steps separated by two bifurcations characterized by the 
coexistence of a continuum of cycles with period 2 and 4 respectively. In Figure |4] 
this corresponds to a raising vertical trajectory in parameter space. As a matter of fact 
this figure contains information on the whole family of bifurcation diagrams like the ones 
in Figures [T] and [2] which stem from, respectively, vertical or horizontal displacements 
in parameter space. Accordingly, a variety of transitions are possible: fixed point to 
chaos, periodic cycle to chaos or fixed point to periodic cycle. We emphasize the role 
the bisecting bifurcations play in mediating these transitions. 

Since piecewise continuous maps are usually analyzed in terms of border collisions 
[3 [S] , it is appropriate to discuss the connection of the bisecting bifurcations with that 
scheme. In the particular case of the Lozi map, the described period— 2 and period— 4 
attractors involving neutrally stable cycles may be seen as border collisions too. To see 
how this comes about we propose to look again at Figure [8] but going in the reverse 
direction. Starting with systems with a = 2 and progressively decreasing its value 
[n = 20, . . . , 16), one finds unstable period— 4 orbits represented by circles. For a = 1.5 
these cycles collide with the borders, namely the discontinuities of LW(x,?/) (dashed 
lines), originating the two segments made of an infinite set of neutrally stable cycles of 
period— 4. After the collision, further decrease of a gives then rise to stable period— 2 
attractors for a < 1 + b = 1.5. Following the way down, the 2-cycle collides, for 
a = 1 — b = 0.5, with the discontinuities of L^^l (x, y) and we see a sudden and ephemeral 
segment of period— 2 neutrally stable cycles. Finally, for a < 0.5 we find the stable 
fixed points Pi. Notice that in these collisions the location of the borders depends on 
a, except the one at x = which is fixed. 

As regards the chaotic regime, two examples of attractors appear in Figure [9} 
intended to be complementary of Figure [8] The upper panel is for a = 1.55 and shows 
an incipient zigzagged structure in phase space as the result of the destruction of the 
regular attractor a = 1.5 in Figure |8| In the bottom panel, with a = 1.7, the strange 
attractor is apparent. It corresponds to the figure in the original reference This 
l^— shaped structure seems to be the prototype of chaotic attractor with \b\ < 1. Points 
in the dotted pattern originate unbounded orbits. 

So far in this section we have considered \b\ < 1. As has already been mentioned, 
systems with |6| = 1 are area-preserving and their attractors offer a large variety of 
structures. Figure [3] is an instance of this. The system with a = l,b = —1, contains 
doubly neutrally stable period— 3 orbits embedded in chaos. In this regard. Figure 



10 illustrates that this situation corresponds also to a border collision. It shows the 



evolution with a of six different fixed points of L^'^^(x,y) (curves coded as dashed and 
dot-dashed lines). The remaining lines in the plot stand for borders. The collisions 
between fixed points and borders at a = 1 are clearly observed. 

The phase space shown in Figure |3] exhibits a mixed structure of periodic an chaotic 
trajectories clusters. The complexity of the phase space is general for \b\ = 1, even for 
non critical cases as a = 1. Figure 11 shows the phase space, at the same scale as Figure 
[3} of the neighbouring system a = 0.9, b = —1. Figure [T2] shows the intricate phase space 
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structure in a larger scale. It is worth noticing that for conservative systems = 1) 
trajectories fill up densely any window in phase space. Thus, to obtain these kind of 
plots, one has to chose not only a large enough transient but also an adequate number 
of initial conditions and iterates to be plotted. Otherwise, the plot will hardly reveal 
inner structures. 



6. Final comments 



In this last section we mention some works related to the Lozi map which have some 
bearing with our results. 

An interesting study of the dynamics of a Henon-Lozi-type map is carried out in 
[13j . The authors introduce a smooth map which depends on a parameter e. In 
the limit e — )■ their map becomes the Lozi map, and they study which properties are 
preserved after the limiting procedure. In particular, they point out the interesting result 
that, whereas the smooth Lozi-like map presents a period-doubling route to chaos, the 
genuine Lozi map does not. This observation buttresses the interpretation advanced 
in the previous section. Instead of the common infinite sequence of successive period 
duplications, here we can observe a transition either from period 1 to 2 for 6 > 1, or 
none for b < 1, before entering the chaotic regime with the distinctive feature that at 
the bifurcation points there is an infinite set of neutrally stable cycles filling a segment, 
or a sector, in the phase plane. 

The theory of border-collision bifurcations provides a classification and a description 
of the bifurcations caused by the collision of fixed points with boundaries in piecewise 
maps. This classification is based on the linearized version of the map around the 
collision point, usually called the border- collision normal form map [8]. However, when 
the normal form is written for the border collisions present in the Lozi map some of the 
non-degeneracy conditions assumed in the classification theorems do not apply [9]. In 
particular, the requirements related to the non-singularity of the coefficient matrices of 
the normal form map are not fulfilled. We refer to ^ for further mathematical details. 

Important explicit results concerning systems with stable periodic orbits embedded 
in chaos have been reported in the literature. In particular, two works [HI [T3] study in 
depth the area-preserving parameter-free maps 

^n+l = 1 ~ 1/n i \Xn\, Un+l = ^n-, (14) 

The former corresponds to the plus sign above and the latter to the minus sign, which 
is introduced because some calculations are considerably shorter. These two maps are, 
except for a refiection transformation [y — i- — y), Lozi maps with a = ±1,6 = 1. In 
an infinite set of period— 3 cycles, as the lowest period, are found embedded in a 



sea of chaos. They are essentially the ones we have reported in section 3.3 In [T3] 



the analogue set is for orbits of period— 6. In both cases further infinite sets of higher 
periods are studied too. 

In this paper we have revisited the Lozi map with the purpose of showing that 
it presents what we have called bisecting bifurcations: those which are mediated by 
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an infinite set of neutrally stable periodic orbits. We have determined explicitly the 
location of some of them as well as their cycle elements. We have also provided an 
explanation for their existence in both algebraic and geometric terms. We deem that 
particular analyses, like the one we have carried out for the Lozi map, may help to 
enlarge our knowledge about the dynamics of piecewise continuous maps. 
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Figure 1. Bifurcation diagram for b = 0.1. Notice that the bisecting bifurcations 
take place at a = 1 ± 6. 
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Figure 2. Bifurcation diagram for a = 1.5. Notice that the bisecting bifurcations 
take place at 6 = ±(a — 1). 
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Figure 3. Attractor of the system a = 1, b = —1. Black and gray regions stand for 
chaotic trajectories. The white areas contain only periodic orbits with neutral stability. 
In particular, the innermost triangle with right angle vertex at (0, 0) contains only 
period— 3 trajectories. 
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Figure 4. (Color online) Behaviour of the Lozi map in parameter space. White areas 
correspond to unbounded trajectories. The bisecting bifurcations take place at the red 
straight lines, a = 1 ± 6. This diagram does not exhaust all the possible attractors the 
Lozi map may exhibit. 
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Figure 6. (Color online) Example of two-dimensional cobweb diagram for the Lozi 
map (a — 0.6, b = 0.3). The point p G is projected to the sm-faces f{x,y) and 
g{x, y) and these projections are brought back down to M? by means of reflections in 
the bisectrix f[x, 0) = a; and (7(0, y) — y in order to find the new point q. For the sake 
of clarity only one small region of each surface is shown. 
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Figure 7. (Color online) Cobweb diagram for the 2nd iterate of the Lozi map with 
b — 0.1 and a = 1 — b. Only the f^^\x,y) surface is shown. The cobweb trajectories 
(blue) of two different initial conditions {xq = 0.3,0.7 and yo = &(1 ^2^o)/a) have been 
represented. 




Figure 8. (Color online) Superposition of twenty attractors illustrating border 
collisions, b — 0.5 and a = 0.1 rt with n = 1,2, ...,20. Period— 2 and period— 4 
attractors associated to neutrally stable orbits (a — 0.5 and 1.5) correspond to the 
segments (solid pattern). In between, stable period— 2 limit cycles. The dashed lines 
stand for the borders at the bifurcation values (as a matter of fact, the border a; = 
is invariant). The circles stand for a period— 4 unstable orbit plotted for the values 
a = 0.1 n, n = 16, . . . , 20 which collides with the borders at the bifurcation. 
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Figure 9. Two strange attractors of the Lozi map for the values a and b specified 
in the panels. Scales are the same as in Figure [8] what allows to appreciate how the 
regular attractor with a = 1.5 in Figure [8] becomes strange as a increases. Points in 
the dotted patterns give rise to unbounded orbits. 
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Figure 10. (Color online) Period— 3 border collision in the Lozi map, illustrated as 
a sequence of snapshots in phase space. The solid line cc = stands for a border of the 
map, which is independent of a. The dotted lines represent further borders of the map 
ata = l,5=— 1. For other values of a 7^ 1 these lines wander across the plot. The 
dashed and dot-dashed lines stand for the evolution of two different orbits of period— 3 
as the value of a varies in the interval [1,3]. The border collision takes place at the 
bifurcation values a = —b = 1. 
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Figure 11. Phase space of the Lozi map for the values a = 0.9 and 6 = — 1. Notice 
that the scales are the same as in Figure [3) 




Figure 12. Phase space of the Lozi map for the values a — 0.9 and b = 



-1. 



